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ABSTRACT 

We have used the Submillimeter Array to image, at ~1.5" resolution, C 2 H 
A = 3 —)■ 2 emission from the circumstellar disk orbiting the nearby {D = 54 pc), 

~8 Myr-old, ~0.8 Mq classical T Tauri star TW Hya. The SMA imaging reveals 
that the C 2 H emission exhibits a ring-like morphology. Based on a model in 
which the C 2 H column density follows a truncated radial power-law distribution, 
we hud that the inner edge of the ring lies at ~45 AU, and that the ring extends 
to at least ~120 AU. Comparison with previous (single-dish) observations of 
C 2 H A = 4 —)■ 3 emission indicates that the C 2 H molecules are subthermally 
excited and, hence, that the emission arises from the relatively warm (T ~ 40 K), 
tenuous {n « 10^ cm“^) upper atmosphere of the disk. Based on these results 
and comparisons of the SMA C 2 H map with previous submm and scattered-light 
imaging, we propose that the C 2 H emission most likely traces particularly efficient 
photo-destruction of small grains and/or photodesorption and photodissociation 
of hydrocarbons derived from grain ice mantles in the surface layers of the outer 
disk. The presence of a C 2 H ring in the TW Hya disk hence likely serves as 
a marker of dust grain processing and radial and vertical grain size segregation 
within the disk. 
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Introduction 


Models of protoplanetary disks orbiting low-mass, pre-main seqnence (T Tauri) stars 
predict that gas-grain processes and high-energy stellar irradiation are important ingredi¬ 


ents in determining disk chemical composition and dust grain evolution (Walsh et al. 2014 


Dutrey et al. 2014, and refs, therein). These evolutionary processes, in turn, have profound 
implications for planet-building scenarios, since disk dust and gas characteristics — such as 
characteristic dust grain sizes, the presence or absence of grain ice mantles, and gas-phase 
molecular abundances — should regulate the rate of planetesimal buildup and the composi¬ 


tions of bodies ranging from comets to giant planet envelopes (e.g., Dutrey et al. 2014, and 
refs, therein). 

Interferometric imaging of molecule-rich disks orbiting nearby T Tauri star/disk systems 


2014 

— provides key insight into these and other planet formation processes (e.g.. 

Qi et al. 

2004 

juris 

Hughes et al. 

2011 

Rosenfeld et al. 

2012 Andrews et al. 2012 ). Recently, for 


example, imaging with the Atacama Large Millimeter Array (ALMA) established that N 2 H+ 
emission from the TW Hya disk displays a central hole of radius ~30 AU (Qi et al. 2013b). 


On the basis of simple gas-phase chemistry, Qi et al. (2013b) interpreted the ~30 AU outer 
edge of the N 2 H+ depletion region as tracing the CO “snow line,” i.e., the innermost radius 
at which CO freezes out onto dust grains. 

In an unbiased molecular line survey of TW Hya spanning the wavelength range 0.85- 
1.1 mm, carried out with the Atacama Pathhnder Experiment (APEX) 12 m telescope, 
we detected strong C 2 H A^ = 4 —)■ 3 emission ( |Kastner et al. 2014). From analysis of 
the hyperhne splitting of this transition, we determined that the emission is optically thick 
and either arises from very cold disk layers or is subthermally excited. We furthermore 
determined a very large beam-averaged C 2 H column density of ~5 x 10^^ cm“^, suggesting 
that, in specihc regions of the outer disk, the abundance of C 2 H in the gas phase may rival 
or perhaps exceed that of CO. As C 2 H is potentially a product of the photodissociation of 
more complex hydrocarbons and dissociative recombinations of molecular ions with electrons, 
its radio emission may trace the influence of stellar irradiation on molecular, protoplanetary 
disks (Henning et al. 2010 Walsh et al.|2010 ). To ascertain the surface brightness distribution 
of C 2 H in the TW Hya disk, which is viewed nearly pole-on [i ^ 6°, assuming a central star 
mass of ~O. 8 M 0 ; Qi et al. 2004; Andrews et al. 2012 ), we have carried out imaging of C 2 H 


X = 3 —)■ 2 emission of TW Hya with the Submillimeter Array (SMA). 
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2. Observations and Data Reduction 

Observations of 262 GHz C 2 H iV = 3 —)■ 2 emission and adjacent continunm centered 
at the position of TW Hya were obtained on 2013 April 5 and 2014 April 8 in the SMA’s 
compact (C) and extended (E) array configurations, respectively, with projected baselines 
ranging from 6.3 m to 164.8 m. The correlator setup included the four brightest hyperfine 
structure transitions of the C 2 H A^ = 3 —?• 2 rotational transition (which he in the 262.004- 
262.067 GHz region; Table 1). The 33^3 — 22,2 and 33^2 — 22,1 transitions were observed in 
lower sideband correlator 104 MHz chunk SOI, and were blended due to the coarse spectral 
resolution of 1.86 km s“^, while the 84^4 — 83^3 and 84^3 — 83^2 transitions were observed in 
lower sideband correlator 104 MHz chunk S02 with spectral resolutions of 0.46 km s“^ (G) 
and 0.23 km s“^ (E), which resolves the two lines (Fig. bottom panel). The observing 
loops used J1037-295 as the gain calibrator, and the bandpass response was calibrated using 
observations of 3G279. Flux calibration was performed using observations of Titan and 
Gallisto. The derived fluxes of J1037-295 were 0.79 Jy on 2013 April 5 and 0.70 Jy on 2014 
April 8 . All data were phase- and amplitude-calibrated with the MIR software package 
Gontinuum and spectral line maps were then generated and GLEANed using the MIRIAD 
software package. 


3. Results and Analysis 


In Fig. we display continuum (266 GHz) and G 2 H emission line maps for the spectral 
region encompassing the 84 ^ — 23^3 and 84^3 — 23^2 transitions of G 2 H. From the merged 
(G-l-E) data, we measure a 266 GHz continuum flux of 737±3 mJy (averaged over both 
upper and lower sidebands) and integrated intensities of 3.74 and 3.25 Jy km s“^ (±0.06 Jy 
km s“^), respectively, for the 84^4 — 23^3 and 84^3 — 23 2 transitions (with absolute calibration 
uncertainties of ~10%). These total flux measurements likely account for all of the 266 GHz 
continuum flux and G 2 H A = 3 —)■ 2 line emission from the TW Hya disk, given the angular 


scales probed by the shortest (~6 m) SMA baselines at this frequency (e.g., Wilner & Welch 


1994). The G 2 H emission is evidently well resolved in the merged G±E conhguration data 


(top row of Fig. [^. The E conhguration map (lower middle panel of Fig. reveals that 
the G 2 H (844 — 23 ^ 3 ) emission has a ring-like morphology, with brightened limbs (ansae) in 
the E-W direction. In the G 2 H ( 84^3 — 23 ^ 2 ) E conhguration emission-line map (lower right 
panel), only the ansae are apparent. 


^http://www.cfa.harvard.edu/~cqi/mircook.html 
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The inner radius of the ring of C 2 H ( 84^4 — 23 3 ) emission in the TW Hya disk, as imaged 
by the SMA, appears to be larger than that of the ring of N 2 H"^ (4-3) emission as imaged 
by ALMA (Qi et ah 2013b). This comparison — which we further explore below (§3.2, 
4.1) — suggests that despite their morphologically similar surface brightness distributions, 
C 2 H and N 2 H~'' may trace different physical and/or chemical processes in the disk. The 
continuum emission is only marginally resolved in these maps, consistent with previous mm- 
wave interferometric imaging observations that established that the continuum emission from 
the TW Hya disk is quite compact relative to the disk molecular line emission (continuum 
and CO outer radii of ~60 AU and ~200 AU, respectively; e.g., Andrews et al. 2012 ). 


3.1. Excitation of C 2 H 

Our analysis of the C 2 H A^ = 4 —)■ 3 hyperhne emission complex had previously indicated 
that the C 2 H emission either arises from very cold (midplane) disk layers, or that the C 2 H 
molecules are subthermally excited (Kastner et al. 2014). The latter interpretation would 


suggest that the relative abundance of C 2 H peaks high in the disk and, hence, is more 


consistent with the chemical simulations of protoplanetary disks presented by Walsh et al. 


(2010), in which the relative C 2 H abundance peaks at relative heights z/R ~ 0.3. 


To investigate whether and how the ratio of C 2 H A^ = 4—)-3toA^ = 3—?-2 line 
intensities as measured with APEX (2.7 Jy km s“^ for the 349.3375 GHz transition; Kastner 
et al. 2014) and SMA (3.74 Jy km s“^ for the 262.0042 GHz transition; §3), respectively, 
constrains the physical conditions (hence height) of the layers of the disk from which the 
G 2 H emission originates, we constructed a grid of models with the web-based RADEX too0 


(van der Tak et al. 2007). RADEX calculates the emergent line flux under the local velocity 


gradient approximation, employing the same molecular data as used for the morphological 
analysis in §3.2. The resulting model G 2 H (4-3)/(3-2) line ratios were then compared with 
the observed ratio, i.e., (4-3)/(3-2) 0.72. The results, which are illustrated in Fig. 

indicate that the observed (4-3)/(3-2) line ratio can be reproduced in two specihc regimes 
of physical conditions: (1) T ~ 15-25 K and logn > 7, or (2) T > 40 K and logn ~ 5-6. 
The latter conditions correspond to subthermal excitation of G 2 H; specihcally, the critical 
densities for excitation of the 3-2 and 4-3 transitions of G 2 H for a gas kinetic temperature 
Tfc = 30 K are Ucr ~0.6xl0^ cm“^ and ~2xl0^ cm“^ respectively (and the values of Her 
increase for smaller assumed values of T^). 

In our tests of models that can reproduce the ring-like morphology of the G 2 H emission 


^http://www.sron.rug.nl/vdtak/radex/index.shtml 
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(§3.2), we find that models in which the C 2 H is assnmed to lie at relatively deep, dense 
layers of the disk consistently overpredict the observed (4-3)/(3-2) line ratio by a factor ~2 
or more, whereas models in which the C 2 H lies high in the disk snrface layers can reprodnce 
the observed (4-3)/(3-2) line ratio to within ~30%. Fnrthermore, the very low excitation 
temperatnre inferred from the C 2 H iV = 4 —)■ 3 modeling {Tex ~ 6 K; Kastner et ah 2014) is 
inconsistent with the observed (4-3)/(3-2) line ratio, nnder conditions of thermal excitation 
(Fig-il. Hence, it appears that the C 2 H is confined to low-n (snrface) disk layers, where it 
is snbthermally excited. 


3.2. Modeling the C 2 H emission morphology 


To constrain the dimensions of the C 2 H ring in the TW Hya disk, we modeled the SMA 


C 2 H data via the same methodology nsed to analyze the ALMA N 2 H’'' (4-3) data (Qi et ah 


2013b). This methodology utilizes a model for the TW Hya disk (described in detail in 


Qi et ah 2004, 2006, 2013b) in which the non-local thermodynamic equilibrium, accelerated 


2D Monte Carlo radiative transfer code RATRAN (Hogerheijde & van der Tak 2000) serves 
to calculate the molecular excitation. The data for the excitation calculations were drawn 


from Schoier et al. (2005), supplemented with H 2 collisional rate coefficients from Spielhedel 


et al. (2012) and electron impact rates appropriate for C 2 H from A. Faure (2014, private 


comm.). As in the Qi et al. (2013b) treatment of N 2 H+, the C 2 H emission is assumed to be 
produced within a vertical disk layer, of constant molecular abundance, whose boundaries 
are expressed in terms of vertically integrated hydrogen column densities as measured from 
the disk surface. The emitting molecular layer is confined to a ring with inner and outer 
edges Rin and Rout-, respectively, with a column density that follows a power-law distribution, 
i.e., N cx RP. 


As noted earlier (§3.1), comparison of the SMA and APEX measurements indicates 
that the C 2 H molecules are snbthermally excited and, hence, reside in tenuous disk surface 
layers. In applying the foregoing ring-like molecular density distribution to model the TW 
Hya C 2 H emission, we therefore focused on simulations in which the C 2 H is confined to a thin 
surface layer. The specific models described here, which are representative of this family of 
simulations, have upper and lower vertical C 2 H boundaries at H 2 column densities of 5 x 10^® 
cm“^ and 8 x 10^® cm“^, respectively. 


With this vertical distribution hxed, we performed hts to the combined (C-l-E) visibilities 
for the resolved C 2 H hyperhne lines, with the ring inner and outer ring edges Rin, Rout and 
the density power-law index p left as free parameters. We hud that the best-ht model, 
obtained by minimizing fhe weighted residuals in the visibilities, has Rin = 45 AU, 
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Rout = 120 AU, and p = —1.8. Based on a comparison of models with the foregoing valnes 
of Rout and p hxed bnt Rin allowed to vary, we hnd that only models with in the range 
40-50 AU well reprodnce the observations (Fig. |^. The onter radius Rout is considerably 
less well constrained, due to degeneracy with the value of p during the model htting (this 


degeneracy between Rout and p is a well-known aspect of power-law disk models; e.g., Mundy 


et al. 1996). We note that these results for the radial distribution of C 2 H within the disk 


i.e., 


the best-£t values for Rout, Rout, and p) are not sensitive to the assumed vertical placement 
of the layer from which the C 2 H emission arises. 

The distribution of C 2 H in the best-£t surface layer model is shown schematically in 
Fig-i This distribution is qualitatively similar to that obtained in the simulations by [Walsh 


et al. (2010), but lies at lower n (in our model, the C 2 H emitting layer corresponds to a disk 
density n ~ 10^ cm“^). The column density and approximate relative abundance of C 2 H 
(with respect to H 2 ) are Nc 2 H = 8.4 x 10^^ cm“^ and XcjH ~ 10“^, respectively, near the 
inner emission boundary at i? ~ 45 AU and ~ 20 AU. The best-fit value of Nc 2 H is similar 
to, though somewhat smaller than, that estimated from measurements of C 2 H 4-3 with 
APEX (Kastner et al. 2014). The value Nc 2 H = 8.4 x 10^^ cm“^ likely represents a lower 
limit on the C 2 H column density at i? ~ 45 AU, given that the C 2 H emitting layer could lie 
somewhat deeper in the disk without violating the constraints imposed by considerations of 
excitation (§3.1). The inferred value of Xc 2 H suggests that the relative abundance of C 2 H 


may rival that of CO in the gas phase, in this region of the disk (see, e.g., |Qi et al. 2013b). 


The result for Xc 2 H is highly uncertain, however, given that the model values of n and T are 
poorly constrained in the tenuous disk layers at high z/R, and that the precise C 2 H emitting 
layer height and thickness is somewhat arbitrary. 


4. Discussion 


4.1. Radial distribution of C 2 H 


The modeling described in §3.1 (and illustrated in Fig|^ indicates that the inner edge 
of the C 2 H ring imaged by the SMA, 45 AU, lies roughly midway between the inner 

edge of N 2 H''' emission as imaged by ALMA (~30 AU; Qi et al.||2013b ) and the outer radius 
of submm continuum emission (~60 AU, as mapped at 870 ptm. and 806 ptm. with the SMA 


and ALMA, respectively; Andrews et al. 2012 Qi et al. 2013b). The N 2 H’'' emission inner 


cutoff likely traces the radius at which CO freezes out onto grain surfaces (Qi et al. 2013b), 


whereas the submm continuum outer cutoff marks the radius at which the large (radius ~1 
mm) dust grain population drops precipitously. These two markers (at 30 AU and 60 AU) 
are well dehned in the respective interferometric data, and neither disk radius is consistent 
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with the C 2 H inner hole dimensions (Fig. [^. The location of the C 2 H emission inner cutoff 
with respect to the CO snow line — and the persistence of C 2 H emission at disk radii 
well beyond the drop in large-grain surface density — suggests that the C 2 H abundance is 
enhanced in radial zones where small, ice-coated grains are plentiful. Indeed, at radii > 60 
AU, such small grains likely dominate the disk’s dust grain mass, based on the comparison 
between the aforementioned submm imaging and the extended (~200 AU radius) scattered- 
light “halo” seen in optical and near-infrared coronagraphic imaging with the Hubble Space 
Telescope (e.g., Debes et ah 2013, see §4.2.3). 


4.2. C 2 H production mechanisms 

In the following, we consider three potential processes that might enhance the abun¬ 
dance of C 2 H in the surface layers of the TW Hya disk beyond ~45 AU: (i) pure gas-phase 
production of C 2 II, involving significant destruction of CO; (ii) X-ray photodesorption of 
organic ices on grains; and (Hi) photodestruction of very small, C-rich grains and PAHs by 
UV and X-rays. As we describe below, mechanisms (ii) and (Hi) may be acting in concert 
to produce the large inferred values of Nc^h- 


4-2.1. Gas-phase production 


A large number of pure gas-phase processes could be responsible for production of C 2 II 
(see, e.g.. Fig. 7 in Henning et ah [2010 ). Many of these channels are initiated via destruction 
of CO and the consequent generation of ionized C (e.g., Aresu et al.j201 H |Cleeves et al.j2013[ ). 
Near the disk surface — above the layer where CO is photodissociated by stellar UV and 
X-ray photons — C ions and atoms can react with H 2 and CH^, respectively, to form C 2 H^, 
followed by recombination and dissociation to C 2 H. This stellar-photon-driven process should 
produce a uniform (as opposed to ring-like) C 2 H column density distribution over the disk 
surface, however, unless the C 2 H-free region within ~45 AU is shadowed from high-energy 
stellar irradiation by material in the inner disk (e.g., in a puffed-up inner rim at ~4 AU; 


Andrews et al. 2012) 


In regions of the disk shielded from UV and X-rays, CR ionization of He results in the 
generation of free C. If the C 2 H ring were generated via such a (CR-initiated) process, then 
it should be conhned to the disk midplane, where high-energy CRs dominate the ionization 


rates (Cleeves et al. 2013). Indeed, our preliminary tests of pure gas-phase chemical models 


(utilizing a physical structure for the TW Hya disk described in detail in Gorti et al. 2011) 


















do yield a ring-like distribution of enhanced C 2 H, with inner radius ~50 AU, deep within 
the disk. However, this placement of the emitting C 2 H (near the disk midplane) appears 
to be inconsistent with considerations of molecular excitation (§3.1). Furthermore, we hnd 
that these pure gas-phase models underpredict the C 2 H column densities by at least an order 
of magnitude (similar modeling results were obtained by Henning et ah 2010). It therefore 
appears that pure gas-phase processes fail to explain the presence of the observed bright 
C 2 H ring in the TW Hya disk. 


4-2.2. X-ray photodesorption of organic ices 


In disk environments such as TW Hya’s, where primitive organics (e.g., CN, HCN and 
C 2 H 2 ) are relatively abundant (Kastner et ah 2014), it is likely that more complex organic 


molecules form and freeze out onto dust grains. Desorption of such ice-coated grains by 
stellar high-energy radiation — especially X-rays — can then be a rich source of dissociation 
products of the parent organics ( Mendoza et al.| 2013). It is not clear, however, whether an 
X-ray-dominated mechanism could produce C 2 H high in the atmosphere of the outer disk. 


given that the stellar X-rays should penetrate relatively deeply within the disk (Cleeves 


et al. 2013). It is also possible that increasing gas/dust ratios in these more remote disk 


regions can facilitate the dredging of small, ice-coated grains to the disk surface, where they 
are efficiently irradiated by ice-destroying UV. However, this UV-driven photodesorption 
process should result in efficient C 2 H production at radii <45 AU, i.e., at least to the CO 
“snow line” at ~30 AU (see also Walsh et al. 2010). Hence — 

§4.2.1) 


as in the case of pure gas- 
- “shadowing” of the disk 


phase processes that are driven by stellar irradiation 
outer regions (by inner disk material) might be required to explain the fact that the C 2 H 
generated by ice photodesorption is evidently restricted to somewhat larger disk radii (~45 
AU). 


4 . 2 . 3 . Photodestruction of small grains and PAHs 


Dust grain fragmentation/coagulation equilibrium models (e.g., Birnstiel et al. 2011 


Gorti et al. 2015) predict both a faster rate of grain growth in the inner disk, where gas 


densities are high, and more efficient collisional shattering of grains at large disk radii, where 
the gas surface densities are lower. The combination of the relatively sharp cutoff of submm 
emission in the TW Hya disk at ~60 AU ( Andrews et al.||2012 ) and the much larger extent of 
the disk in scattered near-infrared light ( Debes et al.||2013 ) is consistent with such modeling. 
This contrast likely reflects a transition within the disk from efficient dust grain coagulation 
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and settling to more efficient grain fragmentation and levitation. 


The outer disk environment hence may harbor an enhanced abundance of very small 
grains and polycyclic aromatic hydrocarbons (PAHs), particularly in the disk’s upper layers 


see, e.g., Fig. 7 in Gorti et al. 2015). Once exposed to stellar UV and X-rays, these small 


grains and PAHs would undergo photodegradation to form loose C 2 H 2 groups (Allain et al. 


1996), which would then photodissociate to produce the observed C 2 H. Preliminary calcu¬ 


lations, based on a simple scheme whereby UV irradiation breaks PAHs into C 2 H 2 units, 
suggest that to explain the C 2 H column density inferred here (i.e., Nc 2 H ~ 10^^ cm“^), the 
relative abundance of PAHs per H nucleus in the upper layers of the disk would need to 
be perhaps an order of magnitude larger than the canonical relative PAH abundance in the 
interstellar medium (~8 x 10“^; Tielens||2008 ). As noted, such an enhanced PAH abundance 
could be a byproduct of the steeply falling gas surface density in the outer disk, which should 
yield disk conditions favorable for the production of small grains. Furthermore, inclusion 
of X-rays in the PAH photodestruction model would lower the required PAH abundance. 
Hence, in light of the modeling results presented in (§3.1, 3.2), we consider photodestruction 
of small grains and PAHs in disk surface layers to be a plausible source of the observed C 2 H. 
The C 2 H production rates from irradiated small grains would be further enhanced if organic 
ices coat these grains (§4.2.2). 


4.3. Origin of the C 2 H ring 


We propose the following schematic model to explain the C 2 H spatial distribution, and 
its apparently subthermal excitation, within the TW Hya disk. Large surface densities are 
required to form submm-sized grains, which also drift radially as the disk evolves; such grain 
growth and migration manifests itself as a sharp edge in the submm dust grain distribution 


(Andrews et al. 2012). A corresponding increase in small grain surface density in the disk 
upper layers moving radially outward, evidently beginning at radii as small as ~45 AU, then 
provides “fertile ground” for UV and soft X-rays generated at TW Hya (via chromospheric 


activity and accretion shocks; e.g., Kastner et al. 2002) to efficiently destroy small grains. 


PAHs, and/or grain mantles composed of organic ices. These photodestruction and pho¬ 
todesorption processes generate hydrocarbons (predominantly C 2 H 2 ) and, via subsequent 
photodissociations, C 2 H. Although pure gas-phase reactions (catalyzed by high-energy stel¬ 
lar irraditation) could also be responsible for C 2 H production, it is unclear whether such 
processes could produce the observed ring-like C 2 H morphology at the large abundances and 
low densities inferred from the data. 


Additional observations and modeling of C 2 H emission from other, nearby protoplane- 
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tary disks are necessary to test and expand on the mechanisms proposed here to explain the 
spatial distribntion and production of C2H in the TW Hya disk. Indeed, few radio measure¬ 
ments of C2H toward disks, either interferometric or single-dish, have been published thus 


far, and the scant interferometric data that do exist (e.g., Henning et al. 2010) are not of 
sufficient quality to determine disk C2H emission morphologies. The results of our recent 
single-dish line surveys of V4046 Sgr (Kastner et al. 2014) and LkCa 15 ( |Punzi et al. 2015) 
indicate that these two disks are particularly deserving of followup interferometric imaging in 
C2H. In each case, the inferred C2H abundance and excitation is similar to that of TW Hya, 
indicating that both disks may also display ring-like C2H surface brightness distributions. 

We further note that in PAH-rich, UV-irradiated environments such as the Orion Bar, 
C2H is merely the most abundant of a series of gas-phase hydrocarbons that are potentially 


detectable in the submm, including C4H, C-C3H2, and C-C3H, among others (e.g., Cuadrado 


et al. 2015). Of these, only C-C3H2 has thus far been detected in a circumstellar molecular 


disk (HD 163296; Qi et al. 2013a). Sensitive searches for this molecule and other simple 
hydrocarbons are warranted for TW Hya and the other two disks just mentioned, to further 
test the various C2 production scenarios proposed in §4.2. 


5. Conclusions 

The SMA imaging reported here reveals that C2H(3-2) emission from the TW Hya disk 
exhibits a ring-like morphology with inner and outer radii of ~1" and ~3", respectively. 
Radiative transfer modeling indicates that the inner edge of the C2H ring lies at ~45 AU 
and that the ring extends to at least ~120 AU; we estimate a lower limit on the C2H column 
density of Nc 2 H ~ 10^^ cm“^ near the inner edge of the ring. The inferred inner cutoff of the 
C2H emission lies roughly midway between the CO “snow line” at ~30 AU (as inferred from 
N2H+ emission line imaging) and the sharp decline in submm-sized dust grains at ~60 AU 
(as inferred from submm continuum emission). Comparison with single-dish observations of 
C2H(4-3) suggests that the C2H emission emanates from tenuous upper layers of the disk 
(characterized by disk heights z/R ~ 0.5 and gas densities n « 10”^ cm“^), where the 
molecules are subthermally excited. While we cannot rule out the possibility that the C2H 
arises in colder (T ~ 20 K), higher-density regions of the disk, it is more difficult for such a 
model to explain the ensemble of observations obtained to date. 

We propose that the C2H emission most likely traces particularly efficient photo-destruction 
of small grains and/or photodesorption and photodissociation of hydrocarbons in the surface 
layers of the outer disk. The presence of a C2H ring in the TW Hya disk therefore may be 
a byproduct of dust grain processing and migration, which is in turn associated with recent 
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or ongoing giant planet bnilding interior to the ring. Given snch a scenario, it remains to 
explain why the inner cntoff in C2H at ~45 AU evidently lies interior to the outer cutoff 
in the large grain surface density. These results hence stress the need for higher-resolution 
(i.e., ALMA) imaging to conhrm and extend the results presented here for TW Hya, and to 
establish the radial and vertical distributions of C2H within additional examples of nearby, 
evolved protoplanetary disks orbiting actively accreting T Tauri stars (e.g., V4046 Sgr, MP 
Mus, T Cha, LkCa 15; Sacco et ah 2014; Punzi et al. 2015). 


The Submillimeter Array is a joint project between the Smithsonian Astrophysical Ob¬ 
servatory and the Academia Sinica Institute of Astronomy and Astrophysics and is funded by 
the Smithsonian Institution and the Academia Sinica. We gratefully acknowledge Alexander 
Faure for providing electron impact rates for C 2 II, and we thank the anonymous referee for 
helpful comments and suggestions. This research is supported by National Science Founda¬ 
tion grant AST-1108950 to RIT and NASA Origins of Solar Systems grant NNX11AK63 to 
SAO. 
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Table 1: SMA Correlator Setups for TW Hya C2H Observations 

C2H transition 

Freqnency 

(GHz) 

Resolntion 
(km s“^) 

iV = 3 ^ 2, J = 7/2 ^ 5/2, F = 4 ^ 3 (34,4 - Sg.a) 

N = 3^2, J = 7/2^ 5/2, F = 3^2 (34,3 - 33,2) 

N = 3^2, J = 5/2^ 3/2, F = 3^2 (33,3 - 22,2) 

N = 3^2, J = 5/2^ 3/2, F = 2 ^ 1 (33,2 - 22,1) 

262.004266 

262.0064034 

262.0648433 

262.0673312 

0.46 (C), 0.23 (E) 
0.46 (C), 0.23 (E) 
1.86 

1.86 
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Fig. 1.— Upper two rows of panels: spectral line maps of the TW Hya disk in 266 GHz 
continnum (left) and C2H (3-2) hyperfine lines (center and right) synthesized from data 
obtained in merged compact+extended (C+E; top panels) and extended (E; middle panels) 
SMA antenna conhgurations. In each panel, the shaded oval indicates the FWHM and 
orientation of the synthesized beams, and the cross represents the position of TW Hya 
(J2000 coordinates a = 11:01:51.834, 5 = —34:42:17.149). Contours levels are 50%, 75% and 
90% of the peak in each image. For the C+E (E) conhgurations, these peak huxes are: 570 
(401) mJy beam“^ for the 266 GHz continuum map, 1.33 (0.49) Jy km s“^ beam“^ for the 
34^4 — 23^3 transition map, and 1.18 (0.41) Jy km s“^ beam“^ for the 84^3 — 23^2 transition 
map, with noise levels of 5 (8) and 0.04 (0.08) for the continuum and line maps, respectively. 
Bottom panel: spectrum of C2H (3-2) hyperhne lines in the 262.005 GHz spectral region. 
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Fig. 2.— RADEX calculations of C2H (4-3/3-2) ratios (see §3.1) as functions of disk H2 
number density (n) and gas kinetic temperature (T) for an assumed C2H column density of 
4 X 10^^ cm“^ (corresponding to a radial position of ~60 AU, in the surface layer model; §3.2) 
and linewidth 0.7 km s“^. The contours indicate the values of n and T needed to match the 
observed C2H (4-3/3-2) line ratio (0.72) and its approximate range of uncertainty (~ ±0.1). 



Fig. 3.— Comparison of observed (left panel) vs. model (remaining panels) maps of C2H 
emission from the TW Hya disk. In each panel, the shaded oval indicates the FWHM and 
orientation of the synthesized or model beam, and the cross represents the position of the 
central star. The inner radii of the model C2H rings are indicated in each panel. All models 
displayed here have the same (best-fit) values of outer radius (120 AU) and density power-law 
index (—1.8). 














